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We investigate the phenomena of the dynamical chiral symmetry breaking ( DCSB ), color su- 
perconductivity ( CSC ), and Bose-Einstein condensation ( BEC ) in a supersymmetric ( SUSY ) 
vector-like SU (N c ) gauge model at finite temperature and density. Both the N = 1 four-dimensional 
and M = 2 three-dimensional cases are considered. We employ the N = 1 four-dimensional gener- 
alized SUSY Nambu— Jona-Lasinio model ( M = 1 generalized SNJL4 ) with a chemical potential 
as the model Lagrangian. The M = 2 three-dimensional theory is obtained by a simple dimensional 
reduction scheme of the four-dimensional counterpart. In order to realize the DCSB and BCS-type 
■ CSC in this model, we introduce a SUSY soft mass term. After adopting the method of SUSY 

auxiliary fields with the Fierz transformation in color and flavor spaces, we discuss several possible 
breaking schemes of the global symmetries of the model. The integrations of the auxiliary fields 
5-H . of composites in the effective potential are performed by using the steepest descent approximation. 

Under the finite-temperature Matsubara formalism, the gap equations are derived and solved. The 
roles of both the boson and fermion sectors in the BEC, DCSB and CSC are examined by the 
quasiparticle excitation spectra and the gap equations. The physical properties of the DCSB and 
the CSC are studied in detail. Examination on the physical property of the BEC is beyond scope 
of this paper. It is found that the BEC, DCSB and CSC can coexist under a condition of model 
parameters. Several important observations are obtained in the method of the construction of the 
SUSY BCS-type theory, starting from a SUSY field theoretical framework. 
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O I- INTRODUCTION 

Based on the successes of the explanations for experiments of superconductivity, the method of BCS theory ( 
Bardcen-Cooper-Schrieffer ) [1] has widely influenced in theoretical physics until now. Through the work of Nambu and 
Q-i Jona-Lasinio ( NJL ), which established the concept of spontaneous symmetry breaking [2], the BCS-NJL mechanism 
was applied to particles [3,4,5], nuclei and hadrons [6,7], and condensed matter physics. In this framework, particle 
and/or hadron masses arc interpreted as results of dynamical chiral symmetry breaking ( DCSB ). Recently, theory 
_C • of color superconductivity ( CSC ) in quantum chromodynamics ( QCD ) has been extensively studied [6]. In high- 
\ density QCD, quarks at the Fermi surface have very high momenta, and SU (3 C ) color gauge coupling becomes weak, 
the gauge interactions between quarks could be satisfactorily handled by perturbation theory at ep > Aqcd [6] 
( ep ; the Fermi energy, Aqcd', the renormalization-group-invariant infrared energy scale of QCD ): The dominant 
interaction between quarks is provided by a one-gluon exchange. Several QCD-motivated NJL-type models are used 
as an effective theory of QCD to investigate hadron physics, meson and diquark mass spectra, and CSC until now [6,7]. 
The theory of CSC motivates us to have interests on several cold dense matters, with promoting developments of 
the methods of field theory at finite temperature and density. In such kind of problems, it might be useful for us to 
employ both the methods of group-theoretical classifications of possible states ( orderings ) , and microscopic theories 
of dynamics of elementary fields of systems. Microscopic theories are parametrized by several parameters, and they 
might reflect a hierarchical structure of Nature. 

Let us assume the existence of supersymmetry ( SUSY ) in the universe. Consequently, all of known elementary 
particles have their superpartners. Phenomenologically, if the Af = 1 SUSY exists, it is a clear fact that it is 
spontaneously broken in the universe. The modern particle theory supposes that the superpartners might have their 
masses in the TeV energy region [8], and now it is a target of the LHC experiment. A natural question arises that, 
what will happen when a condensed matter has a very high-density and its chemical potential or the Fermi energy ( 
the characteristic energy scale of a system ) reach at the TeV energy scale. In such a problem, both particle theory 
and condensed matter theory should be employed. The purpose of this paper is to examine physics, especially the 
DCSB and (color) superconductivity in a SUSY condensed matter with the Fermi energy comparable with masses of 
superpartners. 



Recently, the investigations of CSC from the viewpoint of the SUSY gauge theories were given in literature [9,10,11]. 
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Reference [9] discussed the symmetry breaking patterns of the Af — 1 SQCD ( supersymmetric QCD ) with a nonzero 
chemical potential and a SUSY breaking mass, according as the results of the nonpcrturbativc method of SUSY gauge 
theories [12]. An extension of Ref. [9] to Af = 2 case was given in Ref. [11]. The authors of these works intended 
to classify the possible breaking patterns of the gauge and flavor symmetries of QCD in CSC states, by utilizing the 
solutions of SQCD. However, any gauge-symmetry-breaking two-body pairs like diquarks were not treated in their 
models. The authors of Ref. [10] proposed a toy model for giving diquark-likc condensates. They introduced an 
SO(N) gauge interaction between quark superfields stronger than that of SU(3 C ), and then they argued that the 
SO(N) gauge dynamics gives two-body pair condensations in their model. This SO(N) gauge interaction seems 
artificial. Compared with these works, our motivation and/or context of this paper are different. In this paper, we 
will propose a method for constructing a BCS-NJL-type theory of SU(N C ) DCSB and CSC. We intend to make our 
method parallel with the BCS-NJL theory. We will study DCSB and CSC in a SUSY condensed matter system. 
For this purpose, an SNJL ( supersymmetric Nambu— Jona-Lasinio )-likc model Lagrangian suitable to describe a 
DCSB and a CSC is introduced. We consider the situation where we can neglect the SU(N C ) gauge dynamics, and 
an attractive interaction is provided from an SNJL-type four-body interaction. Historically, the SNJL model was first 
introduced to investigate the phenomenon of DCSB in SUSY field theory [13,14], and an SU(2) x [/(l)-gauged SNJL 
was applied to the top quark condensation theory of electroweak symmetry breaking in the minimal SUSY Standard 
Model [4]. An SU(N C ) SNJL model was used to describe phenomena of phase transitions in the early universe [15]. 

This paper is organized as follows. In Sec. II, the formalism for our investigation is presented. The SUSY auxiliary 
superfields of the composites are introduced. The effective potential and the gap equations for DCSB and CSC will 
be obtained by the method of steepest descent approximation. We use the component field formalism, because it 
is convenient to introduce the SUSY Nambu notation of superconductivity [16,17,18,19,6], and also suitable for the 
method of the steepest descent for our evaluation of the one-loop effective potential. The quasiparticle excitation 
energy spectra are examined in Sec. III. Because we treat a model with a finite chemical potential, a Bose-Einstein 
condensation ( BEC ) can take place in zero modes of scalar fields [9,16,20]. We will find the condition for the 
realization of BEC in our model, and examine the relation between the BEC, DCSB and CSC. Our examination on 
the BEC is only kinematically ( by looking at quasiparticle excitation energy spectra ) and a dynamical description 
of it is beyond scope of this paper. Section IV is devoted for our study on the DCSB without CSC: We solve 
the gap equation for self-consistent chiral mass. Especially the critical coupling of the DCSB is studied. Both the 
Af = 1 four-dimensional and Af = 2 three-dimensional cases are considered. In Sec. V, we examine a CSC of an 
SU(3 C ) x C/(2/)l x U(2f)B, model without DCSB through numerical calculations of the gap equation. The summary 
of this work, discussions, and issues for further investigations are presented in Sec. VI. 



II. FORMALISM 



In this section, we show the formalism for our investigation. First, we introduce the following Lagrangian of an 
Af = 1 SU (N c )-gauged generalized-SNJL model in four-dimensional spacctimc: 

K = Q f e v Q + &e- v Q, 

K x = G 1 (Qh v Q){Q^- v Q) + G 2 (Qh v Q)(Q<fe v Q) + G 3 {&e- v Q)(&e- v Q). (1) 

Here, W a ( W a ) is a chiral ( antichiral ) gauge field strength. Q ( Q ) is an N c x V/-matrix-valued chiral matter 
superfield belonging to the fundamental ( antifundamental ) representations of both SU(N C ) and SU(Nf). The 
definitions are given in terms of component fields: 

Q = 4>(y) + \/29^{y) + 96F{y), Q = <p(y) + V29tp(y) + 99F(y), = a? + i8a*6. (2) 

The chirality constraint D a Q = D a Q = L» Q Q t = D a & = ( D a = d/d8 a + iv^J^d*, D a = -8/89" - i9 a <7^ a d* 
) is satisfied. N c > 3 and Nf > 2 with 3V C > Nf will be considered here. In this paper, we follow the conventions 
for metric, gamma matrices, and spinor algebra given in the textbook of Wess and Bagger [21]. K n + Ki is a 
Kahlcr potential, a real- valued function of the matter fields. K\ includes four-fermion interactions. V = V a Xa ( 
A = 1, • • • , N% — 1 ) denotes a real vector multiplct. The gauge transformation is defined as follows: 

Q^e- lG Q, Q^Q^e lG \ Q — > e iG Q, &^&e~ iG \ 

e ^e- lGt eV G , e- v ^e lG "e- v e- lG , V - V + i(G - G+), Gee^G^. (3) 

A=l 
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Similarly, the global flavor rotation will be given by 

Nj-l 

Q^e~ iT Q, Qt^Qt e fft ) Q^ e ^Q i Q\^Q^ e -iv\ r = ^ r 7 T 7 . (4) 

1=1 

The Hermitian generators Xa of SU(N C ), and Tj of SU(Nf) ( I = 1, • • • , N'j — 1 ) are in the fundamental represen- 
tations, and satisfy 

tr\ A \ B = k c S AB , trT/Tj = k f Su, ^ (AaAa)»j = ^ (T/T/)^- = ^ (5) 

Here, fc c and fc/ are normalization factors taken as k c = kf = 1/2. For the completeness relations, we should also 
include both Ao = l c /\ / 2N c and To = 1// ^/27V/ ( the normalized unit matrices in color and flavor spaces ). Several 
Abelian symmetries are determined as follows: 

U(1) B : Q^e M "Q, Q - e^Q, £7(1) A : Q — > e iaa Q, Q -> e m "Q. (6) 

t/(l)s is a baryon number symmetry, C/(1)a is an axial symmetry. The theory at the classical level has the following 
symmetries of unitary groups: 

SU(N C ) x SU(N f ) L x SU(N f ) R x £7(1) A x U(1) B , (7) 

( = SU(N C ) x [SU(Nf) x J7(1)]l x [SU(Nf) x ?7(l)]ij ), and the quantum numbers of the matter fields are summarized 
as follows: 

Q ■ (n c , N f , Q : (¥ c ,l,iV7,l,-l). (8) 
A chemical potential \i is introduced by 

V -» F + 2^cr 6»® A ®T . (9) 

Namely, /it acts as a zeroth-componcnt of vector associated with the conserved charge of U(1)b [9]. With this form, 
all of the scalar </>,<j> and spinor ip,ip fields share the same /i. <?§ is a bare gauge coupling constant defined at a 
cutoff energy scale. G\, G2 and G3 are coupling constants, and they have a mass dimension [Mass] -2 . Our model 
is non-renormalizable due to K\. In this model, we have introduced a holomorphic soft SUSY breaking mass m~ g for 
gauginos and a non- holomorphic universal soft mass A for squarks [22]. Obviously, these soft masses do not break 
the symmetries given in (7). An i?-symmetry 8 — ► e iar 6 is explicitly broken by including the gaugino soft mass. 

When rrig = A = G\ = G2 = G3 — 0, our model becomes SQCD, while QCD will be obtained by nig — > 00, A — > 00, 
G\ = G2 = G3 = 0. By complete neglection of gauge fields, model (1) becomes a generalized-SNJL ( the original 
version of the SNJL model will be obtained by G2 = G3 = ), and NJL model will be derived by taking A — > 00 ( 
A ^> A ) in the SNJL. Therefore, our SU( A c )-gauged SNJL model (1) might have a rich phase structure. The vacua 
of SQCD were globally determined [12]: SQCD has D-flat directions and the gauge group is broken partially or totally 
along these directions. The origin of the classical moduli space where all of expectation values of squarks vanish, is a 
singular point and the gauge group is totally symmetric. In this paper, we consider the following situation: Starting 
from the origin of the classical moduli space of SQCD, we turn on the interaction of K\. Our model has different 
Kahler metric and Killing potential from that of SQCD. However, the classical moduli space of our model is still the 
same with SQCD because they share the same condition for D-flat directions. Now, we take both the soft masses 
and chemical potential finite, and then we consider a quantum theory. In the case when the one-loop beta function 
coefficient 3iV c — Nf is positive, the SU(N C ) gauge interaction becomes sufficiently weak from both the asymptotic 
freedom in a high-momentum region and a screening effect of a finite-density medium, compared with the interactions 
arising from K x . We should also mention that the Cooper instability of the Fermi surface occurs with any attractive 
interaction. The quantum dynamics of the system is dominated by the attractive interaction of K\, and we restrict 
ourselves on the examination of it. This situation might be achieved by a sufficiently large chemical potential /i with 
the following hierarchy, 

TOg>A>^~0(A) > A (S )qcd, (10) 

where, A is a cutoff for the interactions in K\. A(s)QCD should be determined in SQCD with softly broken M = 1 
SUSY [22]. Under the condition (10), any intermediate processes of the gaugino field do not contribute in the dynamics 
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of our theory. Then the infrared phenomena of the SQCD itself, namely instanton effect, gaugino condensation, etc., 
will be irrelevant. Hence, we do not consider the nonpcrturbative generation of the Afflcck-Dinc-Seibcrg supcrpotcntial 
at Nf < N c in our theory [12]. We skip the question on the origin of the interaction of K\. For example, we could 
imagine that the attractive interaction of K\ arised from the one-gluon exchange interaction in SQCD with the 
condition (10) ( with the complete neglection of gaugino contributions ), though a serious investigation for it demands 
us another work. Because of the nonrcnormalization theorem, the chiral symmetry cannot be broken dynamically if 
the theory maintains the J\f = 1 SUSY exactly in the ordinary SNJL model [13,14]. We speculate it is the case in 
our Af = 1 model, and also in the N — 2 three-dimensional theory obtained by a dimensional reduction. After we 
introduced A, the realization of the DCSB and CSC through the BCS-NJL mechanism in this model might be possible 
through the dynamics governed by K\ . Modern particle phenomenology considers A > TeV, and it is interesting for 
us to imagine the situation discussed in the introduction of this paper. The first term of K i might give a dynamically 
generated Dirac mass, while a left-handed and a right-handed Majorana masses might be generated by the second 
and third terms, respectively. These terms include the factor e ±v to show the SU (N c ) gauge invariance of the theory 
obviously. The special case Gi = G3 ( the left-right symmetric case ) will be taken to keep the parity symmetry 
in our model Lagrangian. In fact, the spin-singlet ( Lorentz-scalar ) BCS pairing gap is given by a parity-invariant 
combination of a left-handed and a right-handed Majorana masses [18,19]. 

After dropping the gauge fields, we apply the Fierz transformation to both the color and flavor spaces: 

JV c 2 -lJV|-l 

{Q ] Q)(QQ ] ) = TTTT E E (Q^aT^QXaTjQ), (11) 

K ^f A=0 1=0 
iV c 2 -lJV|-l 

(g t Q)(QQ t ) = rV E E (Q^BrjQ^(Q\ B TjQ), (12) 

KcK f B=0 J=0 

N 2 c -lN 2 f -l 

(&Q)(Q&) = -j-j^ E E (Q^cT K Q^(Q\cT K Q). (13) 

For example, the numbers of irreducible representations in the color space are estimated by N c ® N c = 1 © N 2 — 1, and 
N c (g)N c = \N C (N C - 1) \N C {N C + 1). Here, lN c (N c - 1) and \N C (N C + 1) give dimensions of antisymmetric and 
symmetric representations, respectively. Consequently, our model Lagrangian will take the following form through 
the method of SUSY auxiliary fields [4,14,15]: 

C = ( d 2 9d 2 9 



(1 - A 2 9 2 9 2 ){Q^e 2 ^° e Q + Qh~ 2 ^° e Q) 
+ J2{^ H l)A! H WAi + S(e)S W Ai(^H {1)AI QXaTjQ) + 5 {9)S\ 1)AI (^Hl )AI - QU^gt)} 



+ £{M^ 2)Bjjff(2)BJ + S (e)S i2)BJ ( k -^H {2)BJ QXbT j Q) + 5(9)S\ 2)bj (^Ih1 2)bj Qt AflTj0 t)} 



+ T,{^ H UckH( 3) ck + S(9)S {3)CK (^H (3)CK - QX C T K Q) + SW^A^^ck - &*c?k&)} ■ 

(14) 

S(i)Ai, 5'(2)_bj an d S^ck are SUSY Lagrange multipliers to keep the relations: 

G\ ~ G2 G3 

H(i)Ai = -tttQ^aTiQ, H(2)b.j = -tttQ^b^ jQ, ff(3)CJf = 77-7— QXc^kQ- (15) 

In terms of the lowest component of superfields, these relations give 4>h {1)ai — • • • , etc. Hence, 

fc (1)A; , • • • , are given as squark pair condensates. On the other hand, the scalars of S^ai, ■ ■ ■ do not have such com- 
pact expressions they might show their "compositeness" in terms of components of Q and Q. Because the superfields 
Q and Q are bosonic, QXb^jQ and QXc^kQ have to be totally symmetric by transpositions. 

To examine a dynamical breakdown of symmetries in the Lagrangian (14), we have to integrate out the matter 
fields Q and Q to obtain a "low-energy effective action", the SUSY Ginzburg-Landau ( GL ) functional given in 
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terms of the SUSY auxiliary fields, similar to the well-known procedure of the theory of superconductivity. The 
Nambu-Goldstone scalars and spinors can be examined by the SUSY GL functional. At zero-temperature and far 
below the compositeness energy scale, the SUSY composite auxiliary fields in (14) appear as independent dynamical 
degrees of freedom, with having radiatively generated kinetic terms ( precisely, Kahler potentials ) of the Lagrange 
multipliers S^ai, • • • [4]. While, a finite-temperature effect excites quasiparticles over a BCS energy gap, and they 
might feel dissociation of condensation pairs. The method of the nonperturbative holomorphic effective action [12] 
could be applied to a SUSY GL functional. ( In this context, the method of the nonperturbative holomorphic effective 
action seems to be a kind of group-theoretical classification of the GL theory in superconductivity. ) Several possible 
solutions of it might be classified by this method, however, this way will give a complicated situation as we have 
mentioned in the introduction [10]. In this paper, we use the method of steepest descent, a kind of the Hartree-Fock 
mean field theory, for the integrations of the auxiliary fields of composites. In fact, at least in condensed matter physics 
of four-dimensional spacetime, the mean field approximation such as the BCS treatment provides a nice description 
for superconductivity. The success of the BCS method is one of the bases of the investigations of the theory of CSC. 
We consider it is also the case in our theory in four dimensional spacetime. In three-dimensional case, a mean field 
description for a phase transition becomes bad because of the existence of strong quantum fluctuations, and usually 
a Kosterlitz-Thouless type transition occurs [23]. However, a treatment of the BCS mean field level will be a starting 
point of further investigations on superconductivity in three dimensions. Expanding the Lagrangian (14) in terms of 
the component fields, eliminating the auxiliary fields of chiral multiplcts through their Eulcr equations, and keeping 
only the relevant terms in the steepest descent approximation, we get 




-{d v - i^5^)4> ] {d v + in5 v0 )(j) - (d v + inS v0 )ft (d v - i/j,5„ )4> 
-</> f (A 2 + Mf)<t> - f (A 2 + Mf)4> + ^M£4> + ftM^cf) 
~i^pa"{d v - inS u0 )ip - i^a 1 '{d v + i^<5„ )V> 

+i>M[ip + $AfPi> + ipMfip + i>MPi> + i>M[ip + $MPi>. (16) 

4>s {1)AI i 0s (2 )B7 an d 4>s {3)C K denote the scalar components of the Lagrange multiplier multiplets. We also assumed 
that, all of scalars of the SUSY auxiliary fields are independent on spacetime coordinates while all of the spinor 
components of them are zero. The mass matrices are defined as follows: For scalar fields, 

Mf = |^0 S(1)A/ A A T / | 2 +4|^^ (2)BJ A B T J | 2 , 

A,I B,J 

A,I B,J C,K A, I 

(17) 

and for spinor fields, 

<^E^u>- A ^ Mf ^^^ (2)BJ A B T 7 , Mi = ^0 S(3)CK AcTk. (18) 

A, I B,J C,K 

The matrices i\)M^i\) and ipM^ ip have to be totally antisymmetric under the exchange of spinor, color and flavor 
indices because of the Pauli principle. This gives us some restrictions on the choices of the symmetries of the SUSY 
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composite auxiliary fields. Now, one finds the partition function as follows: 
2 = / 2W+DWx>*D*IIx>05 (1) "2^ 



A, I 



B,J 



X cxp 



J A.I 1 B.J 2 C,K 



C,K 
|2 



(3)CK 



G 3 



J A,/ B,J C,K 



x exp 



i [ d ^{-k ckf c£ 1%^ + E + E J %^)} 

J A,I 1 B,J 2 C,K 3 



—Vbec + 2zlnDetf2B — zlnDetfip 

vw = -(n c )tQ s ( Pi , = o)n c = ^t( M 2 + Mf + A 2 )^ + 0V + ^f + a 2 ) 

Here, we have used the following definitions of the fields: 



n 



n c + n p ^°(z), * 



H= (*,* T ). 



(19) 
(20) 



(21) 



In the expression given above, the condensate IF is separated from II for taking into account the possible BEC in scalar 
fields. IF does not depend on spacetime coordinates. IF ^ when the BEC takes place, otherwise IF = 0. Vbec 
has no quartic term of n c , and describes the BEC instability of a vacuum. if takes the form of the four-component 
Dirac bispinor in the spinor space. S and S are the Nambu notations [16,17,18], they have 8 x iV c x Nf components. 
* itself belongs to the fundamental representations of both SU(N C ) and SU(Nf)L+R- The definitions of the matrices 
fls and Qf are given as follows: 



&B{p v ) = 
Q. F (jp v ) = 



(d v + i^5 v0 ){d v + i/iSvo) - - A 2 



Mi 



i0 + 7 > - Mf iip - M* T ±^p -M(C{1 + i l5 ) - M* T C(1 - 
■MfC(l + h 5 ) - MPC(1 - if - 7 0T M + Mf + Mf*±=p- 



Ft 1-175 



Ft, 



(22) 
(23) 



Z7 2 7° is the charge conjugation matrix. If we take only X and T in all of the 



where, 75 = 7777, and C 
expansions of the Fierz transformations (11)-(13), we obtain the equivalent result with the £/(l)-gauge case given in 
Rcf. [16]. The numbers of eigenvalues of and Qf become 4 x iV c x Nf and 8 x N c x Nf, respectively. 

Next, we have to consider several possible breaking schemes for the global SU(N C ) x SU(Nf)L x SU(Nf)n x 
?7(1)a x [/(l)s-symmetries of our theory. In the relativistic theory of superconductivity, usually we take the parity- 
invariant Lorentz-scalar symmetry for the spin singlet BCS pairing state [6,16,18,19]. In the case of QCD, scalar will 
be preferred than pseudoscalar due to the f/(l) J 4-breaking instanton effect. Our model is not the same with QCD, 
though we will choose the gap as an antisymmetric function in spinor space with keeping the parity invariance. In 
order to take the symmetries of both the Dirac mass and the BCS pairing gap in Qp as scalar, we choose 



MP = Mf 



MP = -Mi 



(24) 



The gap function of MP 7^ — Mf corresponds to a linear combination of the scalar and pseudoscalar Cooper pairings. 



We will drop the pseudoscalar by the constraint Mf ' = — Mf . In the theory of SU(3 C ) CSC, the symmetry of pairing 
gap of quarks will be determined from the consideration on the attractive channel of one-gluon exchange and the 
Pauli principle [6]. Because we choose superconducting order parameters as antisymmetric ( i.e., a Lorentz scalar ) 
in spinor indices, there are two possible cases: (a) symmetric in both color and flavor indices, (b) antisymmetric in 
both color and flavor indices. The attractive channel of SU(3 C ) is given in the antisymmetric 3 representation of the 
decomposition 3 <g> 3 = 3 © 6. For example, in the case of SU(3 C ) x SU{2f)L x SU{2f)n, the symmetry of a Cooper 
pair is expressed as [6] 



(25) 
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Here, q is a Dirac bispinor of quark field. The quark pair is completely antisymmetric with respect to spinor, color 
( a, b, c ) and flavor ( i,j ) indices. The index c is taken arbitrarily in color space. The color symmetry is broken as 
SU(3 C ) — ► SU(2 C ) x [7(1)a 8 - At iV/ = 3, an order parameter of so-called color-flavor locked phase, where color and 
flavor degrees of freedom are coupled together, is energetically favorable: 

N c N f 

(?Ji(C75)<76j) ~ XIX! c u e iji e «-bJ- (26) 

1=1 J=l 

Here, c/j is determined variationally. Such a color-flavor locked phase are generally favorable at Nf > 3 in SU (3 C ) 
QCD [24]. 

To derive the gap equations, we consider the breaking scheme of the special case SU(3 C ) x SU(2j)l x SU(2f)n, 
the most established breaking scheme in theory of SU(3 C ) CSC. Because the chiral group SU{2f)i, x SU{2f)n is 
unbroken by the CSC itself, we can treat the DCSB and the CSC separately, and it greatly reduces the effort for our 
numerical calculations of gap equations. We take only 

Aa = Ao, As = Ac = A3, Ti = T , Tb = Tc = (27) 

( the index c in Eq. (25) is taken to the third direction in color space ) in Eq.(19), and include one if>f/(3 c )-singlet 
meson. Now, the Nambu notation S has 8 x 3 x 2 = 48 components. Hereafter, we use (f>i = 4>s {1)Q0 an d 4>i = §s,i) 3Q ■ 
The boson and fermion determinants in Eq. (19) will be evaluated in the following forms: 

detfi B = [(po-Ef(p))(po + E^(p))(po-E B (p))(po+E B (p))] 2Nf 

x [(po - £ B (p))(p +^(p))(po - £ B (P))(PV + £ B (p))] {N ^ 2)Nf , (28) 

dctn F = [(p -E F (p))(p + E F (p))(p Q -E F (p))(p Q +E F (p))] 4Nf 

x [(po - £ F (p))(p +£ F (p))(p £ F (p))(p +£ F (p))] 2{Nc - 2)Nf . (29) 

All of the eigenvalues of Qp have two-fold degeneracies in spinor space, and these degeneracies relate to the time- 
reversal invariance of the BCS gap function [18]. The bosonic and fermionic quasiparticle energy spectra are evaluated 
to be, 

Ef(p) = v /p2 + | (/)1 |2 + 4|^ 2 |2 + A 2_ M) (30) 
E B (p) = v / P 2 + |^i| 2 +4|^ 2 | 2 + A 2 + A1 , (31) 

£ B (p) = Vp 2 + I^i| 2 + A 2 - M , (32) 
£ B {p) = A / p 2 + |^ 1 |2 + A 2 +Mj (33) 



££(p) = V(v / P 2 + I^il 2 -M) 2 +4|^2| 2 , (34) 



E F {p) = J(VP 2 + \4>i\ 2 +v) 2 +MH 2 , (35) 



S+(P) = vV+RiP-M, (36) 



£ F (p) = V^T^F + m- (37) 

From the forms of these spectra, we confirm \<j>i\ is a dynamically generated Dirac mass, while 2\(j>2\ corresponds to a 
BCS gap function. The quasiparticles of the branches £±, £ F do not participate the BCS-type superconductivity. 

As we have mentioned above, we will consider the case G 2 = G 3 . The effective action under the steepest descent 
approximation becomes 

S e ff = [ d 4 x\-V B EC-^-^-\ + 2zlnDetf! B -zlnDctf! F , (38) 
J L 4Cti 2G2 J 

2Nf 2N f 

\vbec = (E B ( P =o)rj2m 2 + (E B (P = o)) 2 j2m 2 



2 



(N c -2)N f (N c -2)N f 

+ (£?(p = 0)) 2 Yl \<t> C k\ 2 + (£ B (P = 0)) 2 E |^| 2 , (39) 
fe=l 1=1 

n c -(^, ^ ^f. (40) 
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Here, Il c has been divided according to the breaking scheme in the color space. From the form of Vbec, we recognize 
the following fact: When one of the excitation energy spectra of scalar fields has a zero at p — 0, the BEC takes 
place and at least one of components of n c obtains a finite vacuum expectation value. Now, we introduce the 
finite-temperature Matsubara formalism [20] obtained by the following substitutions in our theory: 



2iri 



y 1 - 



Pa -» iu%,iu}„, 



(41) 



where [3 = 1/ksT ( ks\ the Boltzmann constant, T; temperature ). We take fcs = 1 throughout this paper, and 
are boson and fermion discrete frequencies, respectively. Their definitions are u>^ = 2mr/(3 and u)£ = (2n + l)n/fl 
( n = 0, ±1, ±2, • • • ). The effective potential ( thermodynamic potential ) becomes 



Veff(<t>l,<p2) = Vbec + 



10- 12 



ir , \<h? 



4Gi 

4Gi 2G 2 J (2tt) 3 



+ 



4-/3 7 (2tt)3 F (2tt)3 



2G 2 



+ £ B (p) - - 25f(p)) + (AT C - 2)Nj{£B{p) + £ B {p) £l{p) £ F _{p)) 



-^L M l + e -/5<(P))(l + e -^(P)) + ^£ l n (i _ e -^?(P))(i - e -/J^(P)) 
2(A C - 2)JV,_ ln(i + e _^(p ))(1 + e -^f (p) } + 2(^-2)^ ln(1 _ e -^f(P))(i _ e -^(P)) N 



ANf 



(42) 



To obtain the final expression in Eq. (42), the frequency summations were performed. A three-dimensional momentum 
cutoff A was introduced to regularize the integral. The gap equations are derived in the following forms: 



and, 



= 



dV. 



eff 



= 



dV e 



eff 



d\<h\ 



d\cj> 
(i 



d 3 P 

(2tt)3 



ViF+W E + 2 



• tanh -El + (1 



(coth^ + coth|£ B ) 



t )J-tanh^ F 



+ 4. 



A 2 



/ ^[-^—(tanh^+tanhf^) 

' :(cOth^+COth^ 

V 2 + 2 " 



-^ + — -8^1 



1 



d 3 p 
j2nT 







2' 



1 

w 



tanh -El + tanh 



I- 
2' 



v /p 2 + |</> 1 | 2 +4|<fe| 2 + A 2 



(coth 1^ + coth 



(43) 



(44) 



At first glance, both of these gap equations might have quadratic divergences, and the nature of the divergences might 
be influenced by A. These gap equations correctly give their limiting cases at A — > oo. For example, Eq. (44) gives 
the gap equation of the non-SUSY CSC [6,18,19] at A — > oo. Equations (43) and (44) can have nontrivial solutions 
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at least at G\, G 2 > ( attractive interactions ). The baryon density qb, the conjugate of ii, is found to be 

9Veff . 



QB = -- 



a M + 2Nf J J^f { eJ tanh 2 E+ + W tdnh 2 E - 

+ coth ^E B - coth i^E^j 

+ (N C - 2)N f I ( - tanh + tanh ^£ F + coth - coth ^£ B j . (45) 



(46) 



In the zero-temperature case with |</> 2 | =0 and II C = 0, one has 

QB = 2N c N f J -^0(„ - VP 2 + l^il 2 ) = "cNfj^, 

where, 9(x) is the step function defined as follows: 9{x) = 1 for x > and 6{x) = for x < 0. pi? is the Fermi 
momentum, /x coincides with the Fermi energy y/pp + |<^i| 2 at T = 0, and it is determined by the baryon density of 
fermion. Later, we solve Eq. (44) numerically. We completely neglect the temperature dependence of fj,, and use it 
as an external parameter for our numerical calculations. 



III. QUASIPARTICLE EXCITATION ENERGY SPECTRA 

We will solve the gap equations (43) and (44) to study DCSB and CSC without BEC, in Sec. IV and Sec. V. Before 
that, we examine the quasiparticle excitation energy spectra (30)-(37). Because we consider the case \i ^ 0, we have 
to examine the possibility of BEC ( the condition for BEC ) of the scalar sector in our model. Depending on numerical 
values of the model parameters and gap functions, these excitation branches E B and £ B can have zero-points. In such 

a situation, the logarithmic functions ln(l — e~ /3E + ) and/or ln(l — e~ /3£ + ) in V e ff diverge, and the BEC takes place. 
In this paper, we will not study the physical property of the BEC phase in our model, and only examine the phase 
boundary in the model-parameter space. This can be done by the examination of the quasiparticle energy spectra, 
with neglecting the temperature dependence of [i. 

First, we examine the quasiparticle energy spectra (30)-(37) in several limiting cases. The zero-density case at 
H = 0, the energy spectra becomes 





= E B (p) 


E F + {P) 


= E F (p) 




= £ b (p) 




= £ f (p) 



: VP 2 + l^il 2 +4|<fe| 2 + A 2 , (47) 

■■ y/pt + W+MM 2 , (48) 

vV + |<^i| 2 + A 2 , (49) 

vV + l^il 2 - (50) 

On the other hand, when \(j> 2 \ — 0, one finds the spectra of the case of the DCSB at finite density: 

E B (p) = El (p) = vV + |^i| 2 + A 2 T M, (51) 

El (p) = El (p) = vV + I0i| 2 T M- (52) 
In the case only a CSC gap is generated, \<f>i\ = 0, the spectra become: 

E B (p) = VP 2 +4|^ 2 | 2 + A 2 TAi , (53) 

El(p) = V(lflTM) 2 +4|0 2 | 2 , (54) 

£l{p) = VP 2 TA 2 "TM, (55) 

£±(p) = \P\TH. (56) 

In Eqs. (51), (53) and (55), E B {p) and £+(p) can be negative at /j ^ 0. Because of the positiveness of the Bose 

distribution functions l/(e /3B + — 1) and \/{e f3£ + — 1), and /j corresponds to the Fermi energy of the system at zero 
temperature, 



|0i| <M< Vl^il 2 + A 2 (57) 
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has to be satisfied in the case of Eqs. (51) and (52). By using fj, = \Jpp + |0i| 2 , (57) will be rewritten as the condition 
of the external parameters pf and A: 

0<Pf<&- (58) 
Thus, A is the upperbound for pp in the SUSY theory. Similarly, for Eqs. (53)-(56), 

< V=Pf < & (59) 

has to be satisfied. 

In Fig. 1, we show a typical case of the excitation energy spectra of boson E± (p) and fermion E±(p) quasiparticles 
under the CSC state with a non- vanishing chiral mass |</>i|. E+(p) has a minimum at pp — vV 2 ~~ \4>i\ 2 where the 
excitation energy gap 4|(/> 2 | locates. The appearances of the branches E± (p) and £± (p) are the new phenomenon of 
our SUSY theory compared with the non-SUSY relativistic BCS (color)superconductivity [6,18,19]. All of the spectra 
become parallel with the light cone at \p\ — > oo. In this figure, the Bose branches E+ and £+ have energy gaps at 
\p\ = 0: 



E%(p = 0) = Vl^il 2 + 4|0 2 | 2 + A2- A1 , (60) 

£*(p = 0) = vW+A^-M- (61) 

Note that both E+(p) and (p) depend on p 2 as monotonically incleasing functions, and 

Ef(p) > £f(p) > (62) 

is always satisfied ( £+{p) < is unphysical ). Thus, the BEC never take place in the branch E® at |</> 2 | ^ 0, and 
the condition of BEC will be defined by £+(p = 0) = 0. From this equation, we find 



A = vV-|0i| 2 (63) 
as the condition of the BEC. Under this condition, £+(p) ( = Ef(p) when |0 2 | = ) behaves at \p\ — > as 

S?(p) w ^ ■ (64) 

+ W > 2V(A^ )2 + | 01 |2 

Namely, there is no energy gap of the branch £f under the condition (63). We conclude that, the BEC can take 
place at A = \fjji?- ~ |0iF, and the BEC can coexist with the DCSB and the CSC. Always < pf < A has to be 
satisfied. We should mention that, these results on the relations between the BEC, DCSB and CSC depend on the 
breaking scheme of the gauge and flavor symmetries. The results given above are for the breaking scheme of the 
SU(3 C ) x U(2})l x U{2f)n case discussed in the previous section. 

IV. DYNAMICAL CHIRAL SYMMETRY BREAKING 

In this section, we study DCSB without CSC ( |0i| > 0, \ fo\ = ) and BEC. 

A. M = 1, four-dimensional case 

First, we consider the N = 1 case. At T — and \4>2\ = 0, the gap equation (43) becomes 

2GiJV c iV> f A 2 / 6(p-p F ) 1 



1 = p'dp 



o 



VP 2 + l^il 2 vV + I0i| 2 + A 2 
G i N cN } (, / A2 , , , , 2 .. |2l A +v /A 2 + |0i| 2 



= ^5^(Av^T^F-|^| 2 ln: 



-A0Y 2 + |<H 2 + A 2 + (\<fn | 2 + A 2 ) In 



il 

A+ v /A 2 + |0 1 | 2 + A 2 
Vl^il 2 + A2 



/ 2 . n 12 , I , ,2, PF + VPf + 1^1 1 2 A 

-PFy/Pp + l^i I 2 + l^i I In )L r-^-, J. 



(65) 
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This equation is the same with the U(l) case except the difference of the coefficient of G\ [16]. At A — > 0, the 
right-hand side becomes negative, thus there is no nontrivial solution. With taking the limit A — > oo, Eq. (65) gives 
the gap equation of the DCSB of the non-SUSY NJL model at finite density [7], while at pp — it will give the gap 
equation of the zero-density case of the SNJL model [4,15]. To see the relation of several possible phases in the space 
of parameters GiN c NfA 2 , A/A and Pf/A, we examine of the critical coupling for the DCSB. The determination 
equation for the critical coupling (Gi) cr is obtained as follows: 



7T 2 



{GiNcNfA 2 = — = —7=^ r - (66) 

1 v ^ a 7 A / A A2 

Figure 2 shows (Gi) cr N c NfA 2 as a function of A/A. (Gi) cr N c NfA 2 of the non-SUSY case is obtained by 
lim A / A ^ 00 (Gi) cr A^ c A r /A 2 = ir 2 /(l —p 2 F /A 2 ). The denominator of Eq. (66) includes the finite-density effect on 
(Gi) cr . (Gi)cr of pp 7^ is larger than that of pp = 0. The divergence of (Gi) cr N c NfA 2 at a non-zero value of A/A 
( depends on a numerical value of pp ) indicates the existence of the critical soft mass A® r CSB for the DCSB of this 
model [14]. From (58), \i = pp < A has to be satisfied in Fig. 2: When p F = A, the BEC takes place at the mode of 
p = 0. 



B. AT = 2, three-dimensional case 

In Af = 1 SQED 4 ( four-dimensional supersymmetric quantum electrodynamics ), it was concluded that the dynam- 
ical mass generation never occurs by solving the Schwinger-Dyson equation, while it can take place with introducing a 
SUSY- breaking soft mass [25]. The similar situation happens in Af — 2 SQED 3 ( three-dimensional SQED ) obtained 
by a dimensional reduction of SQED 4 . In QED 3 , there are two-component fermion model and four-component fermion 
model [26]. The Af = 2 SQED 3 obtained by a dimensional reduction of Af = 1 SQED 4 is the SUSY counterpart of 
the four-component-spinor version of QED 3 [27,28,29]. In this paper, we follow the method of Ref. [28], and assume 
all of the quantum fields do not depend on one space-direction. The dimensional reduction is done by the following 
substitution in our gap equations (43) and (44): 

f^, (67) 
2ir r c 

or, equivalently \ J p 2 dp — > j- J pdp. Here, we compactify the third direction ( r c ; compactification scale ). Now, the 
gap equation (43) at T = and |(fe| =0 becomes 

1 = 2G 1 N c N f [ k ^ ( e(p-pp) 1 



irr r 



pdp( 



Vp 2 + \0i\ 2 vV + I0i| 2 + a 2 



2GxN c N 



7T7V 



^(Va 2 + |0x| 2 - y/p F + \h\ 2 - v/A 2 + \4>A 2 + A 2 + Vl^il 2 + A 2 ). (68) 



Here, A is a two-dimensional momentum cutoff. At A — > 0, there is no nontrivial solution similar to the case of the 
Af = 1 four-dimensional theory. The critical coupling will be determined by 



{G^crNcNfA tt/2 



1 . /l I A 2 I A p F 



(69) 



Again, we find that the finite-density effect suppresses the DCSB similar to the Af = 1 four-dimensional case. The 
non-SUSY case is obtained by lim A/A ^ 00 (Gi) cr A r c A / A/r c = (n/2)/(l -p F /A). (Gi) cr N c N f A/r c diverges at A — > 0, 
pp — > 0. Figure 3 shows (Gi) cr N c NfA/r c as a function of A/A. In this figure, pp < A has to be satisfied. The 
dependence of (G\) cr N c NfA/r c on A/A shows qualitatively the same behavior with the Af = 1 four-dimensional case. 



V. COLOR SUPERCONDUCTIVITY AT iV c = 3, N f = 2 

In this section, we restrict ourselves to study a CSC state of the SU(3 C ) x U(2f) L x U(2f) R model without DCSB 
( |0i | = 0, |02 1 > ) and BEC. Usually, the superconducting energy gaps are much smaller than /i in various 
superconductors. Here, we assume it is also the case in our numerical calculation for solving Eq. (44), and choose 
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the model parameters G2, A, /i, and A to satisfy 2 1 02 1 / A* *C 1. The BCS gap function 2|02| is quite sensitive to these 
model parameters, we have to carefully choose the numerical values of them ( a kind of fine-tuning ) to get a physically 
reasonable solution: For example, a too large G2 easily gives a solution larger than /z, and it should be regarded as 
an unphysical one. We also have to take into account pf < A for choosing a numerical value of the soft mass to avoid 
the condition for the realization of the BEC, because our theory cannot handle the BEC state as a stable vacuum. 



A. M = 1, four-dimensional case 

In the N = 1 four-dimensional case, our gap equation (44) is the same with the [/(l)-gauge case [16], except the 
coefficient of the coupling Gi- In Eq. (44), the coefficient 2N f ( coming from the breaking scheme which depends on 
N c and Nf ) is absorbed in G2, and the numerical results are qualitatively equivalent with that of U(l) case: The 
numbers of N c and Nf have no specific meaning in our gap equation of CSC. 

Figure 4(a) shows the gap function 2|0 2 | of the CSC state at T = 0, as a function of the coupling constant G2. 
Equation (44) was solved under the five examples: A/A = 0.45, 0.8, 2, 10, and the non-SUSY case. The non-SUSY 
case, obtained by dropping the contribution coming from the bosonic part in Eq.(44), corresponds to the ordinary 
SU(3 C ) x L/(2/)l x U(2f)n CSC. The gap function 2|02| of all examples show qualitatively similar dependence on G2, 
whether a numerical value of A is larger or smaller than A, and they reflect the nonperturbative dependence on G2. 
In all of the cases, the effect of SUSY suppresses the magnitude of CSC gap functions. When the soft mass becomes 
large, the dependence on G2 becomes strong. The result in Fig. 4(a) also indicates the absence of the critical coupling 
of the CSC similar to the ordinary BCS theory [1] : Any attractive interaction gives a superconducting instability also 
in the SUSY CSC. 

In Fig. 4(b), we show the gap function 2\4>2\ a t T — as a function of A under the energy unit /i = 1. When A takes 
a value closed with /i, the divergence of the gap function becomes moderate compared with other examples because 
of the effect of SUSY: 2|0 2 | diverges almost linearly. When A is large enough, 2|0 2 | will diverge almost quadratically 
under A — > 00. 

Figure 4(c) gives 2|(fo| as a function of temperature T. The gap function continuously vanishes at T — > T c ( T c ; 
the critical temperature ), clearly shows the character of second-order phase transition of the superconductivity. In 
all of the examples shown in this figure, the BCS universal constant 2|0 2 (7 1 = 0)|/T C = 1.76 is satisfied [1]. Both 
2\4>i(T — 0)| and T c depend on G2, A, /x and A, though the ratio of them is independent on these parameters: The 
effect of these model parameters are canceled with each other by taking the ratio. In the non-relativistic [1], the 
relativistic [19], and the U{1) SUSY [16] BCS theories, any spin-singlet Lorentz-scalar BCS gap function satisfies the 
universal constant. We find it is also satisfied in our SU(3 C ) SUSY CSC theory. 

We discuss the thermodynamic property of the SUSY CSC. The definitions of the entropy S and the heat capacity 
G are given by 

In our case, Eq. (42) of V e ff will be used to obtain S and G. In our model, the quasiparticle excitation branches 
will be classified into two classes: (A) E±, E±, and (B) £±, £±. When the BEC does not take place, all of the 
branches E±, E±, £± and £^_ have energy gaps, while only £+ is gapless. The thermodynamic characters of S 
and G are determined by simple sums of the contributions of the ideal gas of quasiparticlcs of these branches. To 
examine the contribution of the branches of (A) in thermodynamics, let us recall the energy spectra shown in Fig. 
1. Because the density of states diverges around the energy gap of E+ [19], thermal excitations of E+ dominate the 
contribution of the class (A) in the integrals of S and G, and other branches of (A) give almost no contribution: 
Namely, the contributions of both the entropy S and the heat capacity G coming from (A) are determined by the 
thermal excitations of the quasiparticles of the branch E+. Because the temperature dependence of 2 102 1 is the same 
with the well-known BCS result, the contributions of thermal excitations of the branches of (A) for S and G are 
suppressed by the factor exp[— 2\<f>2(T)\/T]: They will show the same behavior with the usual BCS theory. On the 
other hand, the low-temperature thermodynamics in the branches of (B) is determined by the gapless excitation £±, 
and it is just the free Fermi gas in our model. Similar to the case of CSC in ordinary QCD, the total thermodynamics 
of the system is given by the superposition of the contribution of the quasiparticles of (A) and (B) . As a consequence, 
we conclude that the thermodynamic character of the SUSY CSC without BEC is the same with the non-SUSY CSC. 
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B. M = 2, three-dimensional case 

We examine the CSC gap equation (44) with the dimensional reduction (67). Figure 5(a) shows the CSC gap 
function 2|(/j 2 | of the Af = 2 three-dimensional case at T = as a function of G 2 /r c . The behavior is almost the same 
with Af = 1 four-dimensional case: No critical coupling, nonperturbative dependence on G 2 , and suppression of 2|(/> 2 | 
by a non- vanishing A. 

Figure 5(b) shows 2|(/j 2 | at T = as a function of A, with the energy unit /i = 1. In the case of A/fi = 1.125, the 
divergence of 2|</j 2 | is quite slow ( almost logarithmically ), while 2|(/j 2 | linearly diverges in the A//i = 5 case. The 
non-SUSY case diverges almost quadratically at A — > oo. These differences of the order ( quadratic, linear, logarithmic 
) of the divergence of 2 1 <^>2 1 in several examples are determined by both the dimensions of space and a numerical value 
of A/fi. 

In Fig. 5(c), we show 2|(/j 2 | as a function of T. The result is qualitatively similar with Fig. 4(c): 2|</> 2 | continuously 
vanishes at T — > T c , and the ratio 2|0 2 (7 1 — 0)|/T C = 1.76 is always satisfied. In fact, this ratio is always satisfied in 
the three-dimensional non-SUSY BCS theory [30]. We find it is also the case in the Af = 2 three-dimensional SUSY 
BCS CSC theory. 

VI. CONCLUDING REMARKS 

In summary, we have discussed the BEC, DCSB and CSC in the SU(N C ) x C/(Y/)l x C/(A r /)_R-invariant generalized 
SNJL model at finite temperature and density. The condition for the realization of the BEC has been examined 
especially in the example of the SU(3 C ) x U{2j)l x f/(2/)# case. The gap equations for the DCSB and the CSC have 
been derived, and have been solved for both the Af = 1 four-dimensional and the Af = 2 three-dimensional cases. We 
have found the finite-density effect in the critical coupling of the DCSB in both cases. The effects of the bosonic part 
of the gap equation (44) in the SUSY BCS-type CSC theory have been reviewed in detail. In the CSC, the SUSY 
effects in the gap function 2|0 2 | become significant when A takes a numerically closed value to fi. Similar to the U(l) 
gauge case of Ref. [16], the CSC shows the BCS character even when /j, is close to A: There is no critical coupling 
in the SUSY CSC, and the BCS universal constant 2\(f> 2 {T = 0)|/T C = 1.76 is always satisfied both in the TV = 1 
four-dimensional and the Af = 2 three-dimensional cases. We should emphasized that, these results are nontrivial and 
indicate the universality of the BCS-NJL mechanism, has nothing to do with the numerical value of ^ ( from eV to 
TeV ), dimension of spacctimc ( 3D or 4D ), SUSY or non-SUSY. 

Finally, we would like to make some comments on several issues. It is interesting for us to examine the collective 
modes in the DCSB and the CSC of our SUSY case at finite density. The collective modes or the Nambu-Goldstone 
modes have been studied in the non-SUSY CSC theory, and an anomalous phenomenon of the Nambu-Goldstone 
theorem was found [31]. An interaction between scalars might alter the excitation energy spectra E± and £± of the 
boson sector, as discussed in the Bogoliubov theory of superfluidity [32]. For the examination of this problem, we 
should take into account the quantum fluctuations around a stationary point of V e f / in our model. 

In this paper, we examined the simplest example of the breaking scheme of SU(3 C ) CSC. However, from a phe- 
nomenological point of view, the results cannot be regarded seriously. Our choice of SU(3 C ) gauge symmetry does not 
have a strong meaning, and should be considered as for a convenience of both the examination of the quasiparticlc 
excitation energy spectra and the numerical calculations of the CSC gap equation (44). A more detailed examina- 
tion of the CSC states under several numbers of N c and Nf is an important issue [24]. Moreover, our SNJL-type 
model does not correspond to a naive supersymmetric extension of the QCD-motivated NJL model usually used as 
a low-energy effective theory of QCD. The most important aim of this work is the presentation of a method for an 
examination of the BEC, DCSB and CSC in a SUSY condensed matter of a non-Abelian gauge model at finite density 
and temperature. In the non-SUSY CSC, the SU(3 C ) gauge interaction is the origin of the attractive interaction, and 
it causes the Cooper instability. Therefore, if one uses an NJL- like model to describe the CSC, the model param- 
eters of it should be chosen from the consideration of the running coupling of QCD gauge interaction and hadron 
phenomenology. A coupling constant for CSC should take a numerical value of 0(GeV~ 2 ), while a cutoff will become 
O(GeV) [6,7], and they are determined by a numerical value of the chiral condensate given by (uu), {dd} and the 
decay constant of pion. At the order of this energy scale of the cutoff, there is no SUSY effect in the CSC in the sense 
of the context of this paper. From the energy scale of A ( >TeV ) relevant for particle phenomenology, our theory 
is relatively closer to the electroweak theory, especially the top condensation model [3,4] or the technicolor theory [5] 
than QCD. As a consequence, it is interesting for us to apply of our theory to several technicolor gauge groups. In 
such a case, the situation \x > Atc should be considered in our method. The extensions to several gauge groups of 
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SUSY grand unified theories ( SGUTs ) are also possible [33,34]. 
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FIG. 1: The excitation energy spectra of quasiparticles of the branches E+(p), E^(p), E+(p) and E+(p) under the CSC state 
with non-vanishing chiral mass. We choose \x = 2, |</> 2 | = 0.02, A — 2 and |</>i| = 1. 

FIG. 2: The critical coupling (Gi) cr N c NfA 2 of the DCSB in the AT = 1 four- dimensional case, shown as a function of A/A. 

FIG. 3: The critical coupling (Gi) cr N c NfA/r c of the DCSB in the M = 2 three-dimensional case, shown as a function of A/A. 

FIG. 4: The CSC gap function 2|0 2 | of the JV = 1 four-dimensional 5^(3,,) x £7(2/)^ x t/( 2 /)fl. model, (a) At T = 0, shown 
as a function of G 2 . We set /x = 0.4 and A = 1. (b) At T — 0, given as a function of A. We set \i = 0.4 and G2 = 1. (c) Given 
as a function of T. The model parameters are set as fi — 0.4, A = 1, G2 = 1. 

FIG. 5: The CSC gap function 2|0 2 | of the N = 2 three-dimensional S77(3 C ) x U(2 f ) L x (7(2/)^ model, (a) At T = 0, shown 
as a function of G2/r c . We set n = 0.4 and A = 1. (b) At T = 0, given as a function of A. We set /1 = 0.4 and G 2 /r c = 0.2. 
(c) Given as a function of T. The model parameters are set as fj, = 0.4, A = 1, G2/r c = 0.18. 



Fig.l, CSC-dispersion(c=l,m=2,d=0.02,u=2) 




Fig.2, DCSB-criticalcouplig(4D) 
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Fig.3, DCSB-criticalcoupling(3D) 
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Fig4(a), CSC-coupling(4D,m=0.4,k=l,t=0) 
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Fig4(c), CSC-temperature(4D,m=0.4,k=l,s=l) 
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Fig5(a), CSC-coupling(3D,m=0.4,k=l,t=0) 
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Fig5(c), CSC-temperature(4D,m=0.4,k=l,s=0.18) 
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